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Abstract 

In the present work the generation of zonal flows in collisionless trapped 
electron mode (TEM) turbulence is studied analytically. A reduced 
model for TEM turbulence is utilized based on an advanced fluid 
model for reactive drift waves. An analytical expression for the zonal 
flow growth rate is derived and compared with the linear TEM growth, 
and its scaling with plasma parameters is examined for typical toka- 
mak parameter values. 
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1 Introduction 



The study of the generation and suppression of turbulence and transport in 
tokamak plasmas is still a high priority topic in theoretical and experimental 
magnetic fusion research. In recent studies, the important role played by non- 
linearly self-generated zonal flows for the regulation of turbulent transport 
has been emphasized |T]- [3]. These are radially localized flows (wave- vector 
Q = (q'x^OjO)), propagating mainly in the poloidal direction, which can re- 
duce the radial transport by shearing the eddies of the driving background 
turbulence. 

The turbulence and anomalous transport observed in tokamak plasmas 
is generally attributed to short- wavelength drift-type instabilities, driven by 
gradients in the plasma density, temperature, magnetic fleld, etc. For the 
hot core region of a tokamak plasma, the two main drift-wave candidates 
are the toroidal Ion Temperature Gradient (ITG) Mode and the coUisionless 
Trapped Electron Mode (TEM). 

Accordingly, the generation of zonal flows by non-linear interactions among 
drift waves has recently been studied both analytically [1]- [TT| and numeri- 
cally [12]- [23]. While a substantial effort has been devoted to the study of 
zonal flow generation by ITG modes (see e.g. [5], [7| and [9] for an analyt- 
ical treatment), very little work has been published on zonal flows driven by 
pure TEM [16]- [T7J. In tokamak plasma experiments, TEM are expected 
to play a dominant role in the hot electron regime (Te > Tj), relevant for ex- 
periments with dominant central electron heating [2l]- [25] and in advanced 
confinement regimes with electron transport barriers {rje > rji). The study of 
zonal flows and turbulence driven by TEM is therefore crucial for the assess- 
ment of these regimes. In addition, an estimate of the zonal flow generation 
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close to marginal stability is essential in order to discriminate between the 
TEM and the potentially important Electron Temperature Gradient (ETG-) 
mode [26]- [27] in comparison of experimental profiles against linear thresh- 
olds. The experimental temperature gradients (or inverse scale lengths) are 
typically found to be about a factor of 2 above the linear thresholds, both 
for ITG and TEM. In the Cyclone study [23], this factor was 1.7 for the 
ITG mode. However, the nonlinear upshift, due to zonal flows, increased 
the effective ITG threshold by a factor 1.5, thus bringing it much closer to 
the experimental gradient. Although a comprehensive and quantitative in- 
vestigation of zonal flow generation would require a nonlinear gyrokinetic 
treatment, a qualitative analytical study, based on a reduced set of fluid 
equations, is feasible and also more transparent in terms of physics interpre- 
tation. 

In the present paper, the generation of zonal flows by pure TEM is stud- 
ied analytically in the limit rji = where the ITG mode is suppressed, using 
a reduced fluid model for the trapped electron dynamics. A system of equa- 
tions is derived which describes the coupling between the background TEM 
turbulence, described by a wave-kinetic equation, and the zonal flow modes 
generated by Reynolds stress forces. The qualitative analytical technique 
used here follows closely the WKB analysis employed in [7J and ^ for 
zonal flow generation by ITG turbulence. The purpose of the study is to 
obtain a qualitative estimate of the zonal flow growth rate driven by TEM 
and to compare with an ITG driven case, and in addition examine its scaling 
with plasma parameters for typical tokamak parameter values. 

The paper is organized as follows. In Section II the model equations for 
the ITG/TEM system is presented and a reduced model for TEM turbulence 
is presented. The equations describing the coupling between the background 
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TEM turbulence and the zonal flows are presented in Section III and most 
explicit derivations are put in the Appendix A. In Section IV the results are 
discussed and flnally there is a summary in Section V. 



2 Reduced Model for Trapped Electron Modes 

The description used for the coupled toroidal ITG and collisionless TEM 
system is based on the continuity and temperature equations for the ions 
and the trapped electrons 

' + V ■ {ujVE + njVi^j) + V ■ {ujVpj + rij-yVj) = (1) 
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where n^, Tj are the density and temperature perturbations (j = i and 
i = et represents ions and trapped electrons) and Vj = ve + v^, + vpj + v-^j, 
ve is the E X B velocity, tJ* is the diamagnetic drift velocity, vpj is the 
polarization drift velocity, is the stress tensor drift velocity and qj is the 
heat flux. The derivative is deflned as d/dt = d/dt + PsCgZ x V0 ■ V and 
(f) is the electrostatic potential. In the forthcoming equations r = Tg/Tj, 



t4 = PsCsV/Ln, Ps = Cs/^ce whcrC Cs = ^jTe/nii, ^Ice = cB / TTleC. Wc alsO 

deflne Lf = — {dlnf /dr)~^ , T]j = Ln/ Lt^, uj^jjuj^j = enQj = ^^9jj where R is 
the major radius, = -^it^ and gj represents the variation of ujE,j along the 
fleld line. The geometrical quantities are calculated in the strong ballooning 
limit {6 = 0, Qj = 1). The perturbed fleld variables are normalized as = 
{Ln/ps)e6(f)/Te, n = {1^/ Ps)Sn/no, fj = (L„/p^)5T,/Teo. The perpendicular 
length scale and time are normalized to ps and Ln/cg, respectively. Equations 
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(1) and (2) can now be simplified to 
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Here ft = rtet/no is the fraction of trapped electrons. The Poisson bracket 
is -B] = dA/dxdB/dy — dA/dydB/dx. The system is closed using the 
quasineutrality condition 



Srii — SUe — SUet + SUef- 



(8) 



where a Boltzmann distribution has been assumed for the free electrons. 
After linearizing equations (4)-(7), the dispersion relation for the coupled 
ITG/TEM system is obtained as 
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where 
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Depending on the plasma parameters, the dispersion relation (9) may contain 
0, 1 or 2 unstable modes. For modes propagating in the ion drift direction 
(usually the ITG mode), Ni < Ne, while for modes propagating in the elec- 
tron drift direction (TEM), N^. < Ni. The modes become de-coupled when 
the inequalities are strong. Thus, for N^. » Ni we obtain a pure ITG mode. 
For pure ITG mode physics, the fluid model used here has been found to be 
in good qualitative agreement with a number of gyrokinetic treatments. For 
example, both the r^j-scaling of the ion heat transport [23J and the nonlinear 
upshift of the linear ITG threshold due to zonal flows [23| has been recovered 
by the fluid model [29j. A more comprehensive version of the model, based 
on the full ITG and TE system (Equations 4-7), has been heavily used in 
predictive transport code simulations [30]- [32] of tokamak discharges. The 
simulation results indicate that the model is able to reproduce experimental 
profiles of temperatures and density, inside the edge region, with good ac- 
curacy over a wide range of plasma parameters. In the limit Ni >> N^., we 
obtain a pure TEM. In tokamak plasmas the TEM is expected to dominate 
in the hot electron regime (Tg >> Tj) and/or in regimes with T]e » rji. 
However, for peaked density profiles (small e„), the ion and trapped electron 
responses are strongly coupled {Ni N^) and a density gradient driven TEM 
appears for weak temperature gradients. In the following we will neglect the 
effects of ion perturbations on the TEM and hence only consider electron 
temperature gradient driven TEM (not including the electron temperature 
gradient mode (ETG)). The dispersion relation then takes the form 
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e = (12) 

Equation 11 describes TEMs driven by R/Ltc and suppressed by i?/L„ lead- 
ing to a linear TEM stability threshold in the parameter r}e = Ln/ Ltc- In the 
considered hmit the TEM is fairly symmetrical to the toroidal ITG mode, 
except that effects of finite- Larmor-radius and parallel electron dynamics do 
not appear in the TEM dispersion relation. The solution to Equation 11 is 
given by 

k / 1 n \ 

oor = -f (e(l-en^)-yen^) (13) 

7 = ky\li^ng iVe-Veth) (14) 

where cu — cUr + and the linear stability threshold is given by 

2 f 10 ^Cen9^ e f,r^ 
^e*.= 3-2 + g^e„, + ^ + ^. (15) 

In this regime we can define a reduced model for electron temperature gra- 
dient driven TEM turbulence by retaining equations (4)- (5) for the trapped 
electron fiuid while neglecting ion dynamics (6)-(7) (see Appendix A). The 
effect of the neglected ion dynamics on the linear physics will in the following 
be quantified by comparing the results of the reduced (equation (11)) with 
the complete (equation (9)) dispersion relation. 



3 Zonal Flow Generation 

In describing the large scale plasma flow dynamics it is assumed that there 
is a sufficient spectral gap between the small scale TEM fiuctuations and the 
large scale fiow. The electrostatic potential is represented by 

0(X, X, y, T, t) = T) + 4>{x, y, t) (16) 
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where y, t) is the fluctuating potential varying on the turbulent scales 
X, t and T) is the zonal flow potential varying on the slow scale X, T 
(the zonal flow potential is independent on Y). 

The evolution of the TEM turbulence in the background of the slowly 
varying zonal flow T) can be described by the wave-kinetic equation [5], [33] 
and [34j for the adiabatic invariant A^^ = Ck\(t>k\^ (see Appendix A for a 
derivation of Nk and Ck)- 

d , ^ 9 / ^\dNk(x,y,t) d /-> _^\dNJx,y,t) 

= -fkNk{x,y,t) - ALjNk{x,y,tf (17) 

Here, Vq is the zonal flow part of the ExB drift. In this analysis it is assumed 
that the RHS is approximately zero (stationary turbulence). The role of non- 
linear interactions among the TEM fluctuations (here represented by a non- 
linear frequency shift Au) is to balance linear growth rate, i.e. jk^kix, y, t) — 
Auj Njs{x , y , ty ~ 0. The TEM turbulence is assumed to be adiabatically 
modulated by the slowly growing potential $(X, T). Equation (17) is then 
expanded under the assumption of small deviations from the equilibrium 
spectrum function; A^^ = Nj^ + Nk where Nk evolves at the zonal flow time 
and space scale {fl, Qx, Qy = 0), as 

-t{n- q^Vgx + Ilk) Nk = K-Q^^-Qj^ (18) 
Here Vg^ = duj/dk^ ~ 0, since the effects of electron FLR is neglected. 



The evolution equations for the zonal flow is obtained after averaging the 
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ion-continuity equation over the magnetic flux surface and over fast scales 
and employing quasi-neutrality (equation 8) 

Here we have assumed that the turbulence is dominated by the TEM (rij << 
fiet) and hence only the small scale self interactions among the TEM are 
contributing to the Reynolds stress in the RHS of (20) [35j. Expressing the 
Reynolds stress terms in equation (20) in Nk we obtain 

-m = J d^kh^kyC^^Nk (21) 

The factor Ck defines the relationship between small scale turbulence and the 
wave action density, see Appendix equation A9 for details. Integrating by 
parts in k^^ and assuming a monochromatic wave packet A^^ = Nq5 {k — ko) 
and using equations (19) and (21) gives 

= -qlC.'klNo (22) 

The dispersion relation for zonal flow Q reduces to 



n = iq.ky^Cj;'No. (23) 

Hence, the zonal flow growth rate scales as oc In expressing the zonal 
flow growth in dimensional form making use of equation (23), it is assumed 
that the background turbulence (in the absence of zonal flows) reach the 
mixing length level for temperature gradient driven modes corresponding to 
Te = ^ . We then obtain (see Appendix A) 

n = iq.kyF-^ (24) 
. - ^^' + ''3 (25) 



r]e - Y^^nQ (1 + 
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where is the zonal flow wave number, ky is the TEM wave number, = 
— Y ('C ~ ^nQ^ + fcnfi') and 7^ is the linear TEM growth rate. The function F 
is usually large in regions close to marginal stability due to the denominator 
in equation (25). This is a result of the quasihnear treatment of the TEM 
perturbations (appearing in equation 23) and Tg = l/{kxLTe)- 

4 Results and discussion 

An algebraic equation (24) describing the zonal flow growth rate driven by 
short-wavelength TEM turbulence has been derived. The zonal flow growth 
rates will in the following be calculated and compared to the linear TEM 
growth rates. First, the linear TEM descriptions are compared. In Figure 
1, the solutions to the full system of 4 equations describing the coupled 
ITG/TE system (squares, equation (9)) is compared to the reduced model 
with 2 equations (asterisks, equation (11)) for the pure TEM. The results 
are shown for rji —0, r = 1, e„ = 1.0, kp — 0.3 and ft — 0.5. The rje scalings 
for the TEM eigenvalues are found to be in good qualitative agreement in 
this regime (for fairly flat density profiles); the growth rates are within 20% 
except close to the hnear threshold. The influence of the ion dynamics on the 
TEM stability would be reduced further in the hmit of cold ions (r >> 1). 

Next, the zonal flow growth rate is studied. In Figure 2 the effects of rje 
and e„ on zonal flow growth rate (normalized to the TEM growth rate) are 
displayed. The other parameters are ft — 0.5, rji — 0, t — 1 and k^p — kyp — 
q^p = 0.3. The results for the zonal flow growth rate are shown for e„ = 0.5 
(with r]eth = 1-38, boxes), e„ = 0.7 (with r^et/i = 1.48, rings). There is a 
significant increment in the zonal flow growth rate (normalized to the linear 
TEM growth rate) just above marginal stability. Part of this increment is 
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Figure 1: The solutions for the TEM eigenvalues using the full system of 
4 equations (squares) is compared to the reduced model with 2 equations 
(asterisks). The growth rate and the real frequency with reversed sign (nor- 
malized the electron diamagnetic drift frequency) vs rje is displayed. The 
results are shown for r = 1, e„ = 1.0, r/j = and kp — 0.3. 
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Figure 2: The zonal flow growth rate (normalized to the linear TEM growth 
rate) vs rje with e„ as parameter is displayed. The results are shown for 
e„ = 0.5 (squares), e„ = 0.7 (rings). The other parameters are r]i = 0, t = 1 
and kxP = kyp = q^p = 0.3 and ft = 0.5. 

due to the reduction of the linear TEM growth. In addition, a resonance in 
the analytical expression for the ZF growth rate (see equation 25) is obtained 
just above the linear TEM threshold which further enhances the ZF growth 
in this region. A similar resonance was found for zonal flows driven by ITG 
modes [9j and [TO]. 

Next the zonal flow growth rates (normalized to the TEM growth rate) 
as a function of e„ (= 2L„/Lb) are displayed (derived from equation. 24) 
with ft (the fraction of trapped electrons) as a parameter. In Figure 3, the 
results are shown for ft = 0.5 (asterisks) and ft = 0.7 (boxes). The other 
parameters are as in Figure 2 with 7]^ = 3. Similar to the results in Figure 
2 a strong growth of zonal flow is obtained close to the e„ threshold (for 
{ft = 0.5, enth ~ 0.09) and {ft = 0.7, enth ~ 0.18)). We emphasize the the 
reduced TEM model used here does not include the Vn drive which will be 
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0.5 1 1.5 



Figure 3: The zonal flow growth rate (normalized to the linear TEM growth 
rate) vs e„ (= 2L„/Lb) with fraction of trapped electrons (ft) as parameter. 
The other parameters are as in Figure 2 with rje = 3, The results are shown 
for ft = 0.5 (asterisks) and ft = 0.7 (squares). 

important in the region of small e„ (e„ << 1). 

Figure 4 compares the zonal flow growth rate generated by pure TEM 
turbulence for rje = 3, rji = to that generated by ITG mode turbulence 
for rje = 0, rji = 3, see equation (33) in Ref. [9] (note the difference in the 
definition of r) as a function of r with e„ as parameter. The comparison 
is done assuming equal saturation levels for the ITG mode and TEM cases. 
The results are shown for e„ = 0.5 (diamonds), e„ = 1.0 (plus) and e„ = 1.5 
(squares). In the case of cold ions (small 1/t), the ITG and TEM turbulence 
generates comparable levels of ZF growth, whereas for equal electron and 
ion temperatures, the ITG mode generates significantly larger levels of zonal 
fiows (typically around 2 times larger ZF growth rates are obtained for the 
ITG case). This is due to the nonlinear diamagnetic effects which significantly 
contribute to the ITG driven ZF growth [5]. The relatively weak generation 
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Figure 4: The zonal growth rate generated from TEM turbulence (for 1]^ = 3, 
rji = 0) compared to that generated from ITG mode turbulence (for r]i = 3, 
?7e = 0) as a function of r with e„ as parameter. The results are shown for 
(en = 0.5) (diamonds), (e„ = 1.0) (plus) and (en = 1.5) (squares). 

of zonal flows obtained here by TEM is consistent with recent results from 
nonlinear gyrokinetic simulations of TEM turbulence |17j . 

5 Summary 

The present paper investigates the generation of zonal flows by collisionless 
trapped electron modes (TEM). An algebraic equation which describes the 
zonal flow growth rate in the presence of collisionless TEM turbulence is 
derived and solved numerically in the strong ballooning limit. A reduced 
model for the electron temperature gradient driven TEM is utilized based on 
an advanced fluid model including the trapped electron continuity and the 
electron temperature equations while neglecting the influence of ion pertur- 
bations. The generation of zonal flows is described by the vorticity equation 
and the time evolution of the TEM turbulence in the presence of the slowly 
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growing zonal flow is described by a wave kinetic equation. 

It is found that the reduced TE model (2-equations) is qualitatively able 
to reproduce the linear physics of the full model (4-equations) in the TEM 
dominated regimes where Te >> Tj and/or in regimes where rje >> r]^.- 
For reasonable fiat density profiles, the linear TE growth rates are typically 
within 20% except close to the linear stability threshold. 

There is a significant increment in the zonal flow growth rate (normalized 
to the linear TEM growth) close to the linear TEM threshold where a reso- 
nance in the ZF generation is obtained. This may result in a larger level of 
zonal flow, and consequently a lower level of TEM turbulence, in this region. 

A qualitative comparison of the zonal flow growth rate generated by pure 
TE and ITG mode turbulence shows that the ITG mode generates signifi- 
cantly larger levels of ZF growth except in the case of cold ions (small 1/r) 
where the TEM driven ZF growth is comparable to the ITG case. 

In the present paper the generation of zonal flows has been studied ana- 
lytically by analysing the linear zonal flow growth rates. A complete analysis 
should also involve an assessment of the stability properties of the generated 
flows. This is more suitable for numerical analysis and is outside the scope 
of the present paper. In addition, a more comprehensive comparison of TEM 
and ITG driven zonal flows should include cases where rjf. rji and/or cases 
where the modes are strongly coupled (for peaked density profiles), using the 
full ITG/TEM system. This is left for future work. 

A Adiabatic Invariant in TEM turbulence 

In this Appendix the derivation of the adiabatic invariant in TEM driven 
turbulence is presented. The method has been described in detail in Ref. [5] 



and [S] (and References therein) and only a brief summary is given here. 
From equations 4 and 5 we get, 

(A.l) 
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wave action density A^^ = |^I/fcP we introduce the normal coordinates ^ 



Here, the interaction between the TEM perturbations have been omitted (see 
discussion after equation 16 for this.) It should be noted that the relationship 
between the electrostatic potential and the trapped electron density riet 
is given by = —j^^het (see equation 8). To determine the generalized 

k = 

^etfc + c^kTet, where is to be calculated. Multiplying equation A.2 by ak 
and adding it to equation A.l gives 

(uetk + OikTetk) + ( + (1 + ens' - "fc-^ (^'7e - genfifj + — e„5fafc 

+ (^(^ngak + ft(^n9^ = - T^etk + "fcTeffc 

The normal coordinates are found if the equation is rewritten as in Ref. [9] 
where 

if 2 \ 2 
Vk = + {I + O^ng - akj^[ve- -engj + ^^engak (A.5) 

a, = ^ (A.6) 

which gives 



-2 - (^ngi + s^ng) + tJitng {Tie " Veth) ^ 

= — ^ n — 7^ — \ ^^-^^ 

j[Ve- ^ (1 + e„^j 
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The linear relations between (j)^, n^tk and 
enables one to express and as 

= rietk + Oik^ek ^ 7 —. (Pk (A.9j 

= l^.r = + /t^) l^^l' = Ck\k\' (A.10) 

Afc = -^fe-en^e+^en^?) (A.12) 



2 V 3 



7fc = kyy ^^n9 iVe - Veth) (A.13) 

The equations A.8 - A.12 describe the normal variables \E'fc, the adiabatic 
invariant A^^ and the linear TEM growth rate, respectively. 
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